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Abstract: The production of forward jets of transverse momentum pr{j) ~ Q and 
large momentum fraction Xj et 3> x probes the onset of BFKL dynamics at HERA. 
A full 0(a 2 ) calculation of the inclusive forward jet cross section is presented and 
compared to the expected BFKL cross section. 

Deep-inelastic scattering (DIS) at HERA provides an ideal place to probe strong interaction 
dynamics. One focus of interest has been the small Bjorken-x region, where one would like 
to distinguish BFKL evolution which resums the leading a s lnl/x terms, from the more 
traditional DGLAP evolution equation , which resums leading a s In Q 2 terms. Unfortunately, 
the measurement of F 2 (x,Q 2 ) in the HERA range is probably too inclusive to discriminate 
between the two ||. 

A more sensitive test of BFKL dynamics at small x is expected from deep inelastic scattering 
with a measured forward jet (in the proton direction) and Pr{j) ~ Q 2 [|J. The idea is to 
study DIS events which contain an identified jet of longitudinal momentum fraction Xj et = 
Pz(jet)/ E pro ton which is large compared to Bjorken x. When tagging a forward jet withp r (j) m 
Q this leaves little room for DGLAP evolution while the condition Xj et 3> x leaves BFKL 
evolution active. This leads to an enhancement of the forward jet production cross section 
proportional to (xj e t/x) ap ~ 1 over the DGLAP expectation. 

A conventional fixed order QCD calculation up to 0(a 2 ) does not yet contain any BFKL 
resummation and must be considered a background for its detection; one must search for an 
enhancement in the forward jet production cross section above the expectation for two- and 
three-parton final states. In this contribution we perform a full next-to-leading order (NLO) 
analysis of this "fixed order" background. Such an analysis has become possible with the 
implementation of QCD radiative corrections to dijet production in DIS in a fully flexible 
Monte Carlo program, MEPJET |5J. 

Numerical results below will be presented both for leading order (LO) and NLO simulations. 
The LO 1-jet and 2-jet results employ the LO parton distributions of Gliick, Reya and Vogt |J 
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Table 1: Cross sections for n-jet events in DIS at HERA at order a°, a s , and a 2 . The jet 
multiplicity includes the forward jet which, when required, must satisfy pr(J) > 5 GeV and 
the cuts of Eqs. (j!]]^). The transverse momenta of additional (non-forward) jets must only 
exceed cuts of 4 GeV (first and third column). This requirement is replaced by the condition 
fcf? > 4 GeV in the second column. No cut is imposed in the 1-jet case at O(a s ) and the 
factorization scale is fixed to Q. 

together with the one-loop formula for the strong coupling constant. At 0(a 2 ) all cross sections 
are determined using the NLO GRV parton distribution functions f{xi,iA 2 F ) and the two loop 
formula for a s (fi 2 R ). With this procedure the 2-jet inclusive rate at NLO is simply given as the 
sum of the NLO 2-jet and the LO 3-jet exclusive cross sections. The value of a s is matched at 
the thresholds (1r = m q and the number of flavors is fixed to rif = 5 throughout, i.e. gluons 
are allowed to split into five flavors of massless quarks. 

Unless otherwise stated, both the renormalization and the factorization scales are tied to the 
sum of parton k^s in the Breit frame, (jlr — Hf — \ Ya k^{i) , where (k^(i)) 2 = 2E 2 (l — cos9 ip ). 
Here 9{ P is the angle between the parton and proton directions in the Breit frame. X^^t(0 
constitutes a natural scale for jet production in DIS |7|] because it interpolates between Q, in the 
naive parton model limit, and the sum of jet transverse momenta, when Q becomes negligible. 

We are interested in events with a forward jet with Pt{j) ~ Q and Xj et 3> x and impose 
kinematical cuts which closely model the HI selection ||] of such events. Jets are defined in the 
cone scheme (in the laboratory frame) with AR = 1 and < 3.5. Here rj = — lntan(#/2) 
denotes the pseudorapidity of a jet. Unless noted otherwise, all jets must have transverse 
momenta of at least 4 GeV in both the laboratory and the Breit frames. Events are selected 
which contain a forward jet (denoted "j") in the pseudorapidity range 1.735 < r](j) < 2.9 
(corresponding to 6.3° < 9(j) < 20°) and with transverse momentum p l ^ b (j) > 5 GeV. This jet 
must satisfy 

Xjet = Pz{])/E p > 0.05 , (1) 

0.5 < p 2 T (j)/Q 2 < 4 (2) 

in the laboratory frame. The condition Xj e t ^> x is satisfied by requiring x < 0.004. Additional 
selection cuts are Q 2 > 8 GeV 2 , 0.1 < y < 1, an energy cut of E{1') > 11 GeV on the 
scattered lepton, and a cut on its pseudorapidity of —2.868 < r)(l') < —1.735 (corresponding to 
160° < 9(1') < 173.5°). The energies of the incoming electron and proton are set to 27.5 GeV 
and 820 GeV, respectively. 

Numerical results for the multi-jet cross sections with (or without) a forward jet are shown in 
Table |l|. Without the requirement of a forward jet, the cross sections show the typical decrease 



with increasing jet multiplicity which is expected in a well-behaved QCD calculation. The 3-jet 
cross section in the last column constitutes only about 10% of the 2-jet cross section and both 
rates are sizable. The requirement of a forward jet with large longitudinal momentum fraction 
(xjet > 0.05) and restricted transverse momentum (0.5 < Px(j)/Q 2 < 4) severely restricts the 
available phase space. In particular one finds that the 1-jet cross section vanishes at LO, due 
to the contradicting x < 0.004 and Xj et > 0.05 requirements: this forward jet kinematics is 
impossible for one single massless parton in the final state. 

Suppose now that we had performed a full 0(a 2 ) calculation of the DIS cross section, 
which would contain 3-parton final states at tree level, 1-loop corrections to 2-parton final 
states and 2-loop corrections to 1-parton final states. These 2-loop contributions would vanish 
identically, once x <C Xj e t is imposed. The remaining 2-parton and 3-parton differential cross 
sections, however, and the cancellation of divergences between them, would be the same as 
those entering a calculation of 2-jet inclusive rates. These elements are already implemented 
in the MEPJET program which, therefore, can be used to determine the inclusive forward jet 
cross section, within the cuts discussed above. At 0(a 2 ) this cross section is obtained from the 
cross section for 2-jet inclusive events by integrating over the full phase space of the additional 
jets, without any cuts on their transverse momenta or pseudorapidities. Numerical results are 
shown in the third row of Table |l|. 

The table exhibits some other remarkable features of forward jet events: the NLO 2-jet 
inclusive cross section exceeds the LO 2-jet cross section by more than a factor of four and the 
3-jet rate at 0(a 2 ) is about as large as the 2-jet rate at 0(a s ). The smallness of the LO 2-jet 
compared to the NLO 2-jet inclusive cross section means that at least three final-state partons 
are required to access the relevant part of the phase space. This three-parton cross section, 
however, has only been calculated at tree level and is subject to the typical scale uncertainties 
of a tree level calculation. Thus, even though we have performed a full 0(a 2 ) calculation of the 
forward jet cross section at HERA, including all virtual effects, our calculation effectively only 
gives a LO estimate of this cross section and large corrections may be expected from higher 
order effects. 

The characteristics of forward jet events are demonstrated in Fig. [I] where the transverse 
momentum and the pseudorapidity distributions of the recoil jet with the highest p l ^ b are shown, 
subject only to a nominal requirement of p l ^ b ,p^ > 1 GeV. Here the recoil system is defined 
as the complement of the forward jet, in the final state which arises in the photon-parton 
collision. Almost all forward jet events contain at least one additional jet in the recoil system, 
with p^ b > 4 GeV and, typically, in the central part of the detector. 

In the usual cone scheme final-state collinear singularities are regulated by the AR separa- 
tion cut while infrared singularities and initial state collinear emission are regulated by the pt 
cut. In 7*p collisions the photon virtuality, Q 2 , eliminates any collinear singularities for initial 
state emission in the electron direction and therefore a large kr is as good a criterion to define 
a cluster of hadrons as a jet as its px- The dashed line in Fig. |TJ(a) shows the kr distribution in 
the Breit frame of the recoil jet candidate with the largest k^. Basically all forward jet events in 
this NLO analysis possess a recoil "jet" with k^ > 4 GeV and would thus be classified as 2-jet 
inclusive events in a variant of the cone scheme where the px > 4 GeV condition is replaced by 
a fcy > 4 GeV cut. This observation makes intuitively clear why we are able to calculate the 
1-jet inclusive forward jet cross section with a program designed for the 2-jet inclusive cross 
section at NLO: there exists a jet definition scheme in which all forward jet events contain at 
least one additional hard jet. 
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Figure 1: Characteristics of the highest transverse momentum "jet" in the recoil system, i.e. 
excluding the forward jet. Distributions shown are (a) da/dpx in the lab frame (solid line) 
and da/dkx in the Breit frame (dashed line) and (b) the jets pseudorapidity distribution in the 
laboratory frame. All distributions are calculated at order a 2 . Jet transverse momentum cuts 
have been relaxed to p l ^ b ,p^ > 1 GeV. 



An estimate for higher order corrections may be obtained by comparing to BFKL calcula- 
tions or to existing experimental results. The HI Collaboration has published such a measure- 
ment which was made during the 1993 HERA run with incident electron and proton energies of 
E e = 26.7 GeV and E p = 820 GeV 0. The acceptance cuts used for this measurement differed 
somewhat from the ones described before. Because of the lower luminosity in this early HERA 
run the Xj et cut on the forward jet was lowered to 0.025 and defined in terms of the jet energy 
as opposed to the longitudinal momentum of the jet in the proton direction, 

Xjet = E(j)/E p > 0.025 , (3) 

and the pseudorapidity range of the forward jet was chosen slightly larger, 1.735 < T](j) < 2.949 
(corresponding to 6° < 6{j) < 20°). Scattered electrons were selected with an energy of 
E{1') > 12 GeV and in the pseudorapidity range —2.794 < T](l') < —1.735 (corresponding to 
160° < 9(1') < 173°). Finally the Bjorken-x and Q 2 ranges were chosen as 0.0002 < x < 
0.002 and 5 GeV 2 < Q 2 < 100 GeV 2 . Within these cuts HI has measured cross sections of 
709 ± 42 ± 166 pb for 0.0002 < x < 0.001 and 475 ± 39 ± 110 pb for 0.001 < x < 0.002. These 
two data points, normalized to bin sizes of 0.0002, are shown as diamonds with error bars in 



Fig. |2|. Also included (dashed histogram) is a recent calculation of the BFKL cross section [ID 

As shown before, the MEPJET program allows to calculate the full 1-jet inclusive forward 
jet cross section^ for x Xj et . The LO result is shown as the dash-dotted histogram in Fig. [| 

1 We have checked that also for the kinematical region considered now almost all forward jet events contain 



2 i i i i i , i , i , , , i i , , , , i 

0.0005 0.0010 0.0015 0.0020 



X Bjorken 



Figure 2: Forward jet cross section at HERA as a function of Bjorken x within the HI 
acceptance cuts (see text). The solid (dash-dotted) histogram gives the NLO (LO) MEPJET 
result for the scale choice \j? R = \xp = £(0.5X)^t) 2 with £ = 1. The shaded area shows the 
uncertainty of the NLO prediction, corresponding to a variation of £ between 0.1 and 10. The 
BFKL result of Bartels et al. |10| is shown as the dashed histogram. The two data points with 



error bars correspond to the HI measurement 



and the NLO result is shown as the solid histogram. The shaded area corresponds to a scale 
variation fi R = fi F = £ | fei^r(0j > fr° m £ = 0.1 to £ = 10, and indicates a range of 
"reasonable" expectations for the forward jet cross section at 0(a 2 s ). 



While the BFKL results [10] agree well with the HI data, the fixed-order perturbative 
QCD calculations clearly fall well below the measured cross section, even when accounting 
for variations of the factorization and renormalization scales. The measured cross section is 
a factor of 4 above the NLO expectation. The shape of the NLO prediction, on the other 
hand, is perfectly compatible with the HI results, and not very different from the BFKL curve 
in Fig. 0. At LO a marked shape difference is still observed, which can be traced directly to 



kinematical arguments given in Ref. [11|. Additional details, including a study of the NLO scale 



dependence of the forward jet cross section, can be found there. First NLO studies for forward 



jet production have been presented in Ref. [TI|. For a study of forward jet cross sections with 
the ZEUS detector, see Ref. 0. 

We conclude that the existing HI data show evidence for BFKL dynamics in forward jet 
events via an enhancement in the observed forward jet cross section above NLO expectations. 



at least one second jet with p l ^ b > 4 GeV and k% > 4 GeV 



The variation of the cross section with x, on the other hand, is perfectly compatible with either 
BFKL dynamics or NLO QCD. Since MEPJET provides a full NLO prediction of the 1-jet 
inclusive forward jet cross section for arbitrary cuts and jet definition schemes, more decisive 
shape tests may be possible as additional data become available. 
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